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Core optics components for high precision measurements are made of stable materials, having 
small optical and mechanical dissipation. The natural choice in many cases is glass, in particular 
fused silica. Glass is a solid amorphous state of material that couldn’t become a crystal due to 
high viscosity. However thermodynamically or externally activated stimulated local processes of 
spontaneous crystallization (known as devitrification) are still possible. Being random, these pro¬ 
cesses can produce an additional noise, and influence the performance of such experiments as laser 
gravitational wave detection. 

PACS numbers: 04.80.Nn, 07.60.Ly, 05.40.Ca 


I. INTRODUCTION 

High-precision measurements always face a lot of noises 
and instabilities. The LIGO project mm have to account 
many fundamental sources of fluctuations. Fluctuations 
of temperature, which are translated into displacement 
of the mirror’s surface through thermal expansion (ther¬ 
moelastic noise) [3113] and change of optical path due to 
fluctuations of refraction index (thermorefractive noise) 
|S] combine producing generalized thermo-optical noise 
mm- Better known Brownian fluctuations causing dis¬ 
placement of the mirror’s surface mm and photoelastic 
effect m produced by these fluctuations form the Brow¬ 
nian branch of noises. 

Not all of noise sources are easy to identify. In this 
work we are trying to estimate a noise coming from struc¬ 
tural transformations in material. Fused silica is a glass - 
neither crystal nor liquid. It is one of polymorphic forms 
of silicon dioxide and its internal energy is higher than 
that of crystalline state. The process of glass to solid 
crystallization at temperatures bellow glass transition is 
often called devitrification and was observed during long¬ 
term heating, under high-intensity laser expositionjTT] or 
ballistic impact [12]. It is essential that different states 
have different material parameters (see table [^, specifi¬ 
cally density and refractive index. 

As one can see from the Table |^the density of fused sil¬ 
ica is smaller than that of crystalline quartz. In this way 



Fused silica 

a-quartz 

Stishovite 

Density, g/cm® 

2.20 

2.65 

4.29 

Heat Capacity, J/(kgxK) 

1.052 

0.740 

0.834 

Refractive index 

1.46 

1.54 

1.80 


Table I. Material parameters of polymorphic forms of silicon 
dioxide 
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some contraction of fused silica samples in time should 
be observed if devitrification takes place. This effect 
of contraction is indeed known for glasses and its rate 
for different materials was measured UMl]. From the 
other hand, each event of local crystallization is a dis¬ 
crete event causing small perturbations. The aim of the 
paper is to calculate the influence of possible local crys¬ 
tallization/reordering processes in the bulk of the matter 
on its surface and to calculate the spectral density of the 
surface fluctuations produced by this effect. We also give 
the empirical estimation of the rate at which such pro¬ 
cesses can happen in fused silica suspension fibers and 
use it to And the absolute values of the corresponding 
noise using the approach introduced in m- 

II. NOISE OF COLLAPSING BUBBLES 

We start considering a piece of glass constituting the 
mirror, a small part of which (which for simplicity we 
shall call a bubble) has changed it’s state. This tran¬ 
sition results in a local change of material parameters 
and also in equilibrium state parameters of such a bub¬ 
ble. One of these parameters is obviously the equilib¬ 
rium volume, that is the volume of non-strained matter 
or exactly mbjpc (where mb is the mass of the bubble, 
Pc is the density of the crystal phase). But as a part 
of the bulk matter in the glassy state, this bubble still 
preserves the volume of the previous state’s equilibrium 
volume mb/pg, where is the density of glass. In this way, 
from the difference of the equilibrium densities of the two 
phases we are getting an initial strain. 

We now have a binary system with the first component 
being a crystal bubble with a deformation Ub 

Uh + ^ot rot Ub — grad div Ub = 0, (H-l) 

Ub\t=o=Uo. (II.2) 

and ci^, consequently and are transversal and 
longitudinal speeds of sound in a crystalline or glassy 
state, uq is time-independent initial deformation field. 

The second component is a glass. As a model task 
to understand the influence of such bubbles on the sur- 
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face we are considering half-space in cylindrical coordi¬ 
nates (p, z, 4>). 

u + rot rot u — cf^ grad div u = 0 (II.3) 
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where a is stress tensor, u is mirror deformation, F is the 
initial boundary (equilibrium form of collapsing bubble 
of the first phase), F' is equilibrium boundary of the sec¬ 
ond phase, n_L and n(i_ are unit perpendiculars to those 
boundaries (see fig. [^. In our case we assume the equilib- 



Figure 1. Evolution of a phase transition bubble. The initial 
state bubble with radius a (left). The bubble changes state to 
a phase with equilibrium radius a' and gets strained (center). 
The system transforms to a new, deformed equilibrium state 
(right). 

rium states to be a sphere with the radius a (F S [r = a]) 
for the glass phase and a' = a^pg/pc (F' g [r = a']) 
for the crystalline phase. The initial deformation can 
be obtained from the stationary spherical problem with 
displacement on the boundary: 

rot rot uq — grad div ug = 0- 

Ur\r=a' = a — a ■ (II-5) 


III. QUASI-EMPIRICAL ESTIMATION 

The above described stress-strain problem is still quite 
complex to solve analytically in arbitrary case. However, 
the initial displacement problem can be modeled by a 
thermoelastic problem, introducing thermal fluctuation 
in a bubble in the form 

T{r) = Tie~^^^ (III.l) 


Then the solution of the elastic equation with the heat 
source 

grad div u — — -- rot rot u = a grad T (III.2) 


l + i^' 


2(1+ iy) 


will provide us a good estimate for (II.l )-(II.51, with local 


pressure substituted by a heat source (a is the coefficient 
of thermal expansion, ly is Poisson coefficient). To And 
the parameters Ti and b consider a pure spherical case 
(expansion of a spericall shell). In |17| we can And ex¬ 
pressions for variations of inner and outer spherical radius 
change as shown on Fig. 


Sr 

SR 


11 + izR^ 
31 — ly r'^ 


(III.3) 


Here i? ^ r is assumed. In our case R is the distance of 



the collapsing bubble from the surface of the mirror, SR is 
the surface displacement, r is simply a - the radius of the 
collapsing bubble and Sr is rteq ” the stationary solution 
of the system (II.l l-(H.5l. For a rough estimation it can 
be approximated as 


Sr 




(HI.4) 


From the other hand, (HI.2) has exact solution in this 
spherical case when R^ b: 


1 + iz b"^ 

Sr =a- - 
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(HI.5) 

(HI.6) 


letting us to estimate Ti and b. We then search for a 
solution of (HI.21 in cylindrical coordinates for half-space 


in the form u = ui + grad^ where (/> takes the right part 

of (IhOI) 


Act> = ^aT 
1 — ly 


(HI.7) 


and gives the boundary for ui problem. So the is a 
simple driven Poissonian solution of well known form and 
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Ml problem can be treated as a boundary-driven halfspace 
problem solved in m- The result for the displacement 
field M should be taken on boundary and averaged over 
the profile of a Gaussian beam with radius w to catch 
the measured displacement of the mirror. 

/ 2 ^ 2 j^y 2 

Uz{x,y,0,Xj,yj,Zj) - dS 

TTW^ 

(III.8) 


After some calculations similar to |3] we obtain the aver¬ 
aged surface response on crystallization of a bubble oc- 
curing at coordinates Xj, yj, zj 


6zj{xj,yj,Zj) 


2aTi{l + v)Tr^^^b^x (IIL9) 

J (27r)3 


Note that 2'K^/‘^aTib'^{l + v) « (l - ^/|^) = 

iVa due to ( |nL4| -( |nL^ , eliminating thermodynamical 
parameters. Here we changed 4/37ra^ to Va - the volume 
of a collapsing region, in attempt to generalize formulas 
to an arbitrary bubble geometry. 

The process of noise consists of discreet collapses 


N{t) 

= X! (Hi.io) 

j 


where Xj , yj, Zj , Tj position and time of a collapse and 
H{t) is a Heaviside step function. For the stationary 
Gaussian process we obtain 


{5z{t)) = A 


^5. = A 


6zj{x,y,z)dV / H{t — T)d 

Jo 

X^VawIi{R/w,L/w)t 
Szj{x, y, zYdVH{—uj)H{ui) 
XeV^Is{R/w,L/w) 

WOj'^ 


(HI.ll) 

(HI.12) 


where A is the process rate parameter - number of events 
per second in unit volume, /; and Is are numerical values 
of underlying integrals, R and L - radius and thickness 
of the mirror. The numerically calculated dimensionless 
integrals R and Is are represented on Fig. 

In [TS] a contraction of a silica Fabri-Perot etalon 
was measured. Two mirrors with a diameter of about 
We = 0.66 cm were connected with a Lg = 10 cm long 
tube with outer diameter of i?e = 2 cm. For this ge¬ 
ometry (HI .81 should be changed as the measurement is 
equivalent to averaging over a ring and not a Gaussian 
spot. Furthermore the |1] approach is not precise as it 
uses an assumption of infinite half-space, while here we 
need a long thin cylinder. To overcome this issue we 
use FEM modeling. A cylinder with the above parame¬ 
ters was modeled using Gomsol Multiphysics. Structural 
Mechanics module was used with two different problem 




IV 


w 


Figure 3. Dimensionless integrals Ii for the case of long cylin¬ 
der (analytics that was used in ( |III.13[)) and numerically cal¬ 
culated Is for a mirror (based on (III.9I). 


formulations: a direct boundary load problem and a pre¬ 
scribed temperature problem (thermal expansion node 
under Linear Elastic Material node). The two solutions 
were found identical with respect to force normalization. 
The main results of the modeling are shown on Fig. EE 

From the simulations it follows that the whole region 
of a cylinder can be subdivided in two parts: 1) a sphere 
with the center in the center of a collapsing bubble, 
touching the closest surface of the cylinder and 2) the 
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Figure 5. Averaged 2 -displacement done by a bubble on the 
cylinder axis as a function of rj and Zj for Gaussian averaging 
(top) and surface averaging (bottom). 


Figure 4. Averaged 2 -displacement done by a bubble on the 
cylinder axis as a function of Zj (top) and cylinder radius 
(bottom). The averaging is made with and without Gaussian 
function (red and blue-plus lines). The theory (green-crosses) 
line in the top figure is ( III.9| |. The i?-dependence is taken 
for different bubble depth Zj and is very close to const/f?^ 
dependence (green-cross line) in all cases of surface-averaging. 


rest of the cylinder. The solution inside the sphere is 
close to the spherically symmetric solution (III.3I. Out¬ 
side the sphere the solution is close to constant spherical 
field. The averaged elementary response from one bub¬ 
ble in case of surface averaging is found to be practically 
constant with depth and varies less than by 0.8% with 
offset from the axis. 

The averaged elementary response from one bubble in 
case of Gaussian averaging is in good agreement with the 
theory (IIL9I till the depth of 2.5w, and approaches the 
first modeling case after S.Gic (see fig. |4]top). That can 
be explained easily, having in mind that the cylinder ra¬ 


dius here was R = 3w. The idea is that the bubble “feels” 
only the closest boundary (the one that is touched by the 
earlier mentioned sphere). For the depths less than R 
(and close to the cylinder axis) the governing boundary 
is the front surface, making the problem similar to half¬ 
space and demonstrating appropriate transversal varia¬ 
tions for Gaussian averaging and constant for surface av¬ 
eraging (see fig. [^. For greater depths the governing 
surface is the side surface, and the distance from it to 
the collapses near cylinder axis is constant, providong a 
constant response from depth. 


In this way, for noise calculation we assume that (III.91 
is valid until the depth of the bubble is smaller than R 
and stays constant for larger depths with an error no 
larger than 13%. However, as LIGO mirrors have L k, R 
we thus stay in half-space approximation for spectral 
density and do not need the z > R extension. For the 
time constant determination we assume the elementary 
response from one bubble to be constant and equal to 
6zj{0,0, R). With these arguments the contraction and 
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noise (III.91 can be found analytically. Then the limit of 
small w should be taken to remove the Gaussian beam ra¬ 
dius from expressions. It ensures that the response value 
at the depth R is equal to the one of surface averaging 
case. Thus we get (52j(0,0, i?) oc R~^ which is shown as 
green line on Fig. (right). For the process parameter 
we obtain 


(fe) 2 2Rl 

L, VaiRl-wl 


(III.13) 


We assume that the size of collapse is of the order of 
silica molecule 14 ~ ^^1?^ = 45 x 10“^° m^ (a ~ 0.22 
nm), where Msi 02 = 60 g/mol - molar mass, Na - Avo- 
gadro number. From [15] we get izz = —5.8 x 10“^® per 
second and calculate A « 1.54 x 10^® events per second 
per m^. The resulting spectral density of devitrification 
noise at 100 Hz for ii = L = 20 cm and w = 6 cm is 


Sz — 


z^VaLe Is{R/w,L/w) 


WU}'^ We Il{Re/We,Le/We) 

= 6.31 X 10-2^ m/Hz^/^ (III.14) 


which is 8000 times smaller than the Brownian noise in 
substrates and coatings for LIGO mirrors iHiiinj. 


IV. DEVITRIFICATION NOISE IN STRING 
SUSPENSIONS 


is already 2 x 10^ times smaller then the Brownian sus¬ 
pension noise. 


V. CONCLUSION 


The main uncertainty of our estimates is the average 
radius of collapsing bubbles. Our initial idea to esti- 

AdvLIGO Noise Curve: P. = 125.0 W 
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Figure 6. Creep suspension noises from m using process 
parameter of devitrification (III. 131 and devitrification noise 
(|III.14[) in mirrors together with other Advanced LIGO noises. 


Another estimate of event-based noise was made re¬ 
cently for suspension fibers by Yu. Levin HU. He consid¬ 
ered spontaneous discrete stress relaxation events (creep 
events) in suspension strings. However he also suffered 
from the lack of the process rate parameter and event 
volume values and thus could not obtain absolute noise 
values. One can speculate about the origin and direc¬ 
tion of the creep events, but local reordering may be one 
of the sources. So we can use formulas (49) and (52) 
from m to estimate devitrification noise in suspensions, 
changing R{V'^) (Levin’s notation) to NVgXV^ (our no¬ 
tation), where 14 = ttVsIs ~ the volume of the string and 
N is the number of strings. Three devitrification noises 
together with existing LIGO noises are shown on Fig. 

This naive estimate of suspension noise due to devit¬ 
rification does not take into account complex multistage 
suspension system of LIGO. Furthermore the time pa¬ 
rameter for the loaded case (mirror mass is about 40 
kg) is probably smaller because the extension caused by 
massive mirrors opposes crystallization with contraction. 
Nevertheless we present at least the upper bound, which 


mate it from the crystal-glass internal energy difference 
encountered a serious problem as this energy is also not 
known exactly. Different estimates from literature give 
values varying by two-three times. 

Note that spectral density (HI. 121 strongly depends 
on average radius of collapsing bubbles for constant A: 
Ssz Aa®. However, taking in account that the value of 
\a? is taken from contraction rate ( |HL13| ), we get that 
effectively ^/Ssz i.e. increase of the radius a 

10 times will increase the estimate (HI. 141 30 times. It 
means that reliable knowledge of collapsing radius a (as 
well as A) is very important. 
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